Singularity of Berry Connections Inhibits the Accuracy of Adiabatic Approximation 
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Adiabatic approximation for quantum evolution is investigated quantitatively with addressing its 
dependence on the Berry connections. We find that, in the adiabatic limit, the adiabatic fidelity 
may uniformly converge to unit or diverge manifesting the breakdown of adiabatic approximation, 
depending on the type of the singularity of the Berry connections as the functions of slowly- varying 
parameter R. When the Berry connections have a singularity of 1/R a type with a < 1, the adiabatic 
fidelity converges to unit in a power-law; whereas when the singularity index a is larger than one, 
adiabatic approximation breaks down. Two-level models are used to substantiate our theory. 
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The adiabatic theorem, as a fundamental theorem in 
quantum mechanics, plays a crucial role in our under- 
stand and manipulation of microscopic world[l|. Re- 
cent years witness its growing importance in the con- 
text of quantum control, for example, concerning adia- 
batic transportation of the newly formed matter — Bosc- 
Einstein condensate 0, as well as for adiabatic quantum 
computation^. 

Despite the existence of an extensive literaturejj] on 
proofs of estimates needed to justify the adiabatic ap- 
proximation, doubts have been raised about its validity 
leading to confusion about the precise conditions needed 
to use it [1, 0] . In the present paper, we point out that 
the above confusions can be avoided with formulating the 
quantum adiabatic approximation within parameter do- 
main rather than time domain. Within this formulation, 
we address the fidelity of the adiabatic approximation 
quantitatively (up to the prefactors of power-law or ex- 
ponential forms). In particular, we find that in the adi- 
abatic limit, the behavior of adiabatic fidelity depends 
on the singularity of the Berry connections. That is, 
the Berry connections of the system determine how good 
the adiabatic approximation is. As we know, in adia- 
batic quantum search algorithms, the upper bound of 
adiabatic fidelity (measure of distance between exactly 
solution and adiabatic approximation) is essential to the 
search time 0, H| ■ 

The system we consider is a Hamiltonian containing 
slowly-varying dimensionless parameters R(£) belongs to 
a given regime [Ro,Ri], saying, H(R(t)). Initially we 
have a state, for example, the ground state |i?o(R.(^o))) 
with energy E (R(t )). The wave function \^(t)) fulfills 
the usual Schrodinger equation, i.e., = H (R(t))\I> (t) , 
with h = 1. The above problem has a well-known adia- 
batic approximate solution , 



|* Qd >= e -^ t£ « dt e^|S (R(t)) > 



(1) 



with 70 — i f dt < E \E >, the geometric phase term 
@, E3] • The above equation is the explicit formulation of 
the adiabatic theorem stating that the initial nondegen- 
erate ground state will remain to be the instantaneous 



ground state and evolve only in its phase, given by the 
time integral of the eigenenergy (known as the dynamical 
phase) and a quantity independent of the time duration 
(known as the geometric phase). 

The problems is, how close the above adiabatic ap- 
proximate solution to the real solution |^(t)). To clarify 
the above question and formulate it quantitatively, we 
introduce two physical quantities, namely, adiabatic pa- 
rameter and adiabatic fidelity, respectively. 

The dimensionless adiabatic parameter is defined as 
the ratio between the change rate of the external param- 
eters and the internal characteristic time scales of the 
quantum system (i.e., the Rabi frequency \E m — E n \), 
used to measure how slow the external parameter changes 
with respect to time, 



IRI 



\E n (R) - E m (R)\ 



, 771 71. 



(2) 



e — > corresponds to adiabatic limit. 

Adiabatic fidelity is introduced to measure how close 
the adiabatic solution to real one, F a d = (SE'ft) | VP^d) | 2 . 
The convergence of the the adiabatic fidelity to unit uni- 
formly in the range R G [Ro,Ri] in the adiabatic limit 
(e — » 0), indicates the validity of the adiabatic approxi- 
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Figure 1: Schematic illustration of the quantum adiabatic 
evolution formulated in parameter condition, see text for de- 
tailed description. 
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mation. Evaluating fidelity function will give an estima- 
tion on how good the adiabatic approximation is. 

In Fig.l we schematically illustrate the physical pro- 
cess we describe above. Our main result is that, the 
distance between the adiabatic solution and real one con- 
sists two parts: the fast oscillation term and the secular 
term. The time scale of the oscillation is the Rabi period, 
its amplitude is proportional to the square of the adia- 
batic parameter. The amplitude of the secular term is 
exponential small (~ exp — 1/e) suppose that the Berry 
connections of the system are regular, and turn to be 
power-law (~ e x ,x < 2) if the Berry connections have 
singularity or the external parameters vary in time non- 
linearly. 

We start our statement with writing the wavefunc- 
tion as a superposition of the instantaneous eigenstates, 

!*(*)> = E„c'„(i) e -/ tdt (^-^(«)i^(«)))|i? n (R(t))), 

and suppose initial state is the ground state, i.e., Cb(t = 
0) = 1, C n (t = 0) = 0, (n ^ 0). Then the adiabatic ap- 
proximate solution takes form of Eq.(l) and the adiabatic 
fidelity F ad = \(M!(t)\y ad )\ 2 = \C \ 2 ~ 1 - \AC n \ 2 ,(n ^ 
0) . To evaluate the adiabatic fidelity we need quanti- 
tatively evaluate the change of the coefficients C„ with 
respect to time. 

Substituting the above solution into Schrodinger equa- 
tion, we have following differential equation for the coef- 
ficients, 



dt 



-a nn R) — (E„ 



,R))it 



at 
(3) 

where a nm (ft) is the Berry connection. Both off-diagonal 
and diagonal Berry connections have clear physical mean- 
ing and important applications (llj . We first suppose 
these Berry connections and the gradient of instanta- 
neous eigenenergy are not singular (NS) as the functions 
of the external parameters, i.e., 

cw(R) = (E n \iV n \E m );p n (R) = V R £„(R), NS. (4) 

Under the above condition and in the adiabatic limit, 
to estimate the distance between the adiabatic solution 




and the real one, we need quantitatively evaluate the 
change of the coefficients C n with respect to time. The 
right-hand of Eq.(3) contains unknown C m , to the first 
order of approximation, we take Co = 1 and C m — 0, to ^ 
in the right-hand of Eq.(3). Then, Eq.(3) shows that, 
the change consists two parts: the fast oscillation term 
and the secular term. The time scale of the oscillation 
is the Rabi period, its amplitude is proportional to the 
adiabatic parameter suppose that the Berry connections 
are regular with limitation. Whereas the secular term 
maybe is exponential small of form e"~ or power-law 
depending on the Berry connections in following. 

Let us denote n — j (E n — Eo)dt, and then the upper 
bound of the increment on the coefficients (n ^ 0) can 
be evaluated as follows 



i(H-l) e id n 



,(Ro) E n - E 



Mr. 



a n oRd6 r , 



,(Ro) 



(Hi), 



(5) 
M n (6) 



where we set that in the right-hand of Eq.(3) the coef- 
ficients C m ~ for to ^ 0, and C'o ~ 1 since we want 
to estimate the upper bound of the adiabatic approxi- 
mation. For simple and without losing generalization, in 
following deductions we regards the slowly-varying as a 
scalar quantity, and assume that ^ ~ e 9(R) with func- 
tion g(R) being regular with limitation. 

The main contribution to the first term on the right 
hand comes from the pole point, 9^ — f c (E n — E )dt ~ 
- J °{E n — Eo)/g{R)dR, determined by the equation 
E n (d c n ) — E (8 c n ) = 0. Because the non-degeneracy 
of the system, the solutions of the above pole points are 
complex with imaginary parts. Let ojq be the singular- 
ity nearest the real axis, i.e., the one with the small- 
est (positive) imaginary part (see Fig. 2). Then, the 
first term is approximately bounded by exp(—Imu!o) ~ 

e -X\j\^*c)\ {En _ Eo)/g{R)dR ^ wh . ch . g ^ gecular 

term and is exponential small in adiabatic limit . The 
second integral is proportional to e, which gives the esti- 
mation on the fast oscillation term. 

Now, we consider the situation that Berry connection 
has singularity in the parameter regime [Ro,Ri], e.g., 
R = 0, taking the form We then evaluate the above 
integral in the neighborhood domain [— e, e] of the singu- 
lar point, other regime is regular and can be treated in 
the same way as before. Then, 



IACJ 



jf(E n -E )dt 



a n0 (R)dR\ 



(7) 



Figure 2: The integral path and the singular points in the 
complex plane. 



The situation is divided into two cases: a < 1 and a > 1. 
For a < 1, this type of singularity can be removed be- 
cause the integral is finite. The integral in the neigh- 
borhood domain [— e, e] of the singularity contributes a 
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quantity of order e^ 1_<T \ We thus expect that the adia- 
batic fidelity approaches to unit uniformly in the 2(1 — a) 
power-law of the adiabatic parameter, i.e., 

2(l-<x) 



1 - F ad ~ e z 



(8) 



For the case of a > 1 , the singularity is irremovable and 
the adiabatic approximation is expected to break down. 

The above discussions is readily extended to the case 
that the slowing-varying parameters change nonlinearly 
with time, i.e., R = et a ', with a is any positive number. 
The nonlinear time dependent parameter has many phys- 
ical origins, for example in molecule spin system the ef- 
fective field vary in time nonlinear ly[13j. Another field of 
broad examples is quantum optics, Rabi frequency cou- 
pling different levels (i.e., stimulated Raman adiabatic 



passage) is usually nonlinear dependence on time 14 1. 
Here we suppose that the Berry connections of quantum 
system are regular with limitation as the function of the 
parameter R and the level spacings are of order 1. To 
apply our theory, we introduce the new parameters R' 
and e' through the expressions e' — e° and R 1 = e't. 
Then, R = R' a . As a function of the new parameter R', 
the singularity of the Berry connections determined by 
the 4jg ~ 1/ R' 1 ° . Our discussions are divided into two 
cases: i) a > 1 and ii) a < 1. In the former case, the 
Berry connections as the functions of the new parameter 
are regular, so the adiabatic fidelity is determined by the 
short-term oscillation and is expected to converge to one 
in a power-law of e' 2 . Then we have 



l-F ad 



(9) 



In the latter case, the Berry connections as the functions 
of the new parameter are singular of type 1/R' 1 ° . For- 
tunately, this singularity is removable, it gives an upper 
bound of the adiabatic Fidelity as e' 2rT , i.e., 



1 - F ad 



e 2 . 



(10) 



Notice that in this case, the upper bound of the adiabatic 
fidelity is independence of the nonlinear index a . 

In following, we use two-level models to substantiate 
the above discussions. Our model is denoted by S a : a 
spin-half particle in a rotating magnetic field, its Hamil- 
tonian reads, 

H a {R{t)) = ~ (sin 9 cos f(R(t))a x + sin6smf (R(t))a y 



-cos 9<j z ), 



(11) 



u>o is defined by the strength of the magnetic field, and 
<7j, (i = x,y,z) are Pauli matrices. f(R) is a function of 
slowly-varying parameter R. The above 2x2 matrix is 
readily diagonalized for fixed R and we then obtain its 
instantaneous eigenenergies E± = -t^rr, and instanta- 
neous eigenvectors, 



\E a _(R)) = 



i 2 sin - r 

-e 2 cos ; 



,\E1(R)) 



e 2 cos • 
e 2 sin ^ 



The Berry connections are derived as follows, 



1 a df c 



1 • n d f 
— sin 0——. 

2 dR 



(13) 



i)We first consider following two cases: f(R) — In \R\ 
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Figure 3: Adiabatic fidelity evolves for different type of sin- 
gularity. 
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Figure 4: Upper bound of the adiabatic fidelity in system S c 
for different type of singularity. F£$ n is the minimum fidelity 
in the parameter range R € [— 2n, 2n] 

and f(R) = |i?| 1_<T with u < 1, respectively. And 
here the slowly- varying parameter is supposed to linearly 
change with time, i.e., R(t) = tot. lu is the rotating fre- 
quency of the magnetic field. Apparently, the Berry con- 
nection has singularity of the form 1 / R a at point R = . 
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These systems are complicated and analytic solutions are 
not reachable. We thus make numerical simulations on 
the adiabatic fidelity by directly solving the Schrodinger 
equation with the 4 th — 5 th Runge-Kutta adaptive step 
method. Our results are shown in Fig. 3 and Fig. 4. In 
the Fig. 3, it is clearly shown that, without the singular- 
ity (Fig. 3a), the distance between the adiabatic solution 
and real solution is determined by the fast oscillation, 
therefore gives the upper bound of square adiabatic pa- 
rameter. With the singularity (Fig. 3b), the upper bound 
is determined by the type of the singularity of the Berry 
connections as we discuss above. For the situation that 
the Berry connections have irremovable singularity (i.e., 
(7 = 1, see Fig. 4b) the adiabatic fidelity converges to 0.53 
rather than one implying the failure of adiabatic approx- 
imation. For the case of the removable singulary of the 
Berry connection(cr < 1) the adiabatic fidelity converges 
to unit in the power-law dependence of the adiabatic pa- 
rameters as we expect (see Fig. 4a). 

ii) We then set f(R) = R and R varies in time non- 
lincarly, i.e., R — esign(t)\t\ a . In this case, we see that 
(Fig. 5), for a > 1, the adiabatic fidelity converges to 
unit in a power-law of exponent 2/cr; for tr < 1, it clearly 
demonstrates that the exponent turns to be two, inde- 
pendence of the nonlinear index a. These numerical sim- 
ulations corroborate our theory. 
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Figure 5: Upper bound of the adiabatic fidelity for different 
type of the nonlinearly-varying external parameters. 

In the above discussion we discuss quantum adiabatic 
issue under the parameter conditions. We emphasize that 
the adiabatic problem can only be well formulated in the 
parameter domain but not in time domain. Therefore, 
it requires that the Hamiltonian depends on the time 
only through the slowly- varying parameters (i.e., taking 
form H(R(t))) and, the range of the parameters (i.e., 
(Ro,Ri)) can be reached at certain time (i.e., to and ti) 



for any small adiabatic parameter (e) as sketched in Fig.l. 
Usually, the smaller the adiabatic parameter, the longer 
the time duration (i.e., t\ — to) is needed. In following 
we address this point with an example raised in 0]. It 
is constructed from S a for the case of f(R) = R, R = ujt 
through following relation, 



H ce = _u a \t)H a {t)U a {t) 



(14) 



with U a (t) = Texp(f*H a (t')dt') the time evolution 
operator of system S a . Its explicit analytic expres- 
sion is readily obtained Q. After lengthy deduction, 
we obtain the explicit expression of the Hamilto- 



nian H ce = fL{i) 
2uiu>o cos cos 2 wt/ 2 



a, where L(t) = (sui^Wq + 
f uj 2 cos zut)/zu 2 , ^^-smzut, 



cos 8 



2uuq sin i 



sin (9 sin vat/2), and w 



\/ u>o + oj 2 + 2uiloo cos 9. 

It is easy to verify that L(t) is a unit vector, i.e, 
L(t)| = 1. The eigenvalues and eigenvectors for this sys- 
tem are, 



E c l = ± 



LJO 

2 ' 



H e (L)) = 




where <fi = h arctan(L 2 /£i). 
Berry connections as follows, 



(15) 



We then can obtain the 



2-^3 



2(1 -L 2 Y 2(1 - L\) 1 
L2 L\ 



0), 



(- 



2 v /r~L| 2 v /r~Lf 2 v /r^T| 



.(16) 



As I/3 < 1, the Berry connections are not singular. For 
d|L|/<if ~ w, the adiabatic parameter of this system is 
e = lu/ujo- 

The controversy is that, even though in the adiabatic 
limit e — > 0, the adiabatic fidelity calculated in the time 
domain t 6 [— 27r, 2ir] does not diverge to unitjl]. More- 
over, with changing the sign of the above Hamiltonian 
and re-calculating the adiabatic fidelity in the time do- 
main [— 27r, 2n], we find that the adiabatic fidelity con- 
verges to unit in the adiabatic limit. The above result is 
rather confused. The reason for the above controversy is 
that the problem is discussed in time domain rather in 
parameter domain. 

To resolve the above controversy, we check the above 
system in the parameter domain . First, after trans- 
formation (14), the R — Lot acted as the slowly-varying 
parameter in H a (t) system no longer should be chosen 
as the slowly- varying parameter of the new system H ce , 
because the Hamiltonian H ce depends explicitly on the 
time not only through R. Instead, L(t) can serve as the 
slowly- varying parameters. However, the range of the pa- 
rameters keeps the same order of the adiabatic parameter 
(i.e., |Li(f = t%) — L (f = to)\ ~ to) and tends to zero 
in the adiabatic limit no matter how long the evolution 
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time t\ — io is. This completely counters to our picture 
schematically plotted in Fig.l. Above analysis indicates 
that the system H ce can not be well formulated in the 
parameter domain, it essentially not a system that adia- 
batic theory can applies to. If one discuss the dynamics 
of this system in the time domain as shown the above, 
any strange things can happen. 

In summary, we investigate the fidelity for quantum 
evolution under the parameter domain with addressing 
the adiabatic approximation quantitatively. Within this 
framework, we clarify the confusions in applying quan- 
tum adiabatic theory, and find that the singularity of 
Berry connections inhibit the accuracy of the adiabatic 
approximation. Our estimation on the adiabatic fidelity 
has important meaning in the practical adiabatic quan- 
tum search algorithms . 
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